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Abstract. The validity of a closure called the minimal tau approximation (MTA), is tested in the context of dynamo theory, 
wherein triple correlations are assumed to provide relaxation of the turbulent electromotive force. Under MTA, the alpha effect 
in mean field dynamo theory becomes proportional to a relaxation time scale multiplied by the difference between kinetic and 
current helicities. It is shown that the value of the relaxation time is positive and, in units of the turnover time at the forcing 
wavenumber, it is of the order of unity. It is quenched by the magnetic field - roughly independently of the magnetic Reynolds 
number. However, this independence becomes uncertain at large magnetic Reynolds number. Kinetic and current helicities are 
shown to be dominated by large scale properties of the flow. 
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1. Introduction 

In many branches of astrophysics it is necessary to use turbu- 
lent transport coefficients. This allows one to describe the col- 
lective effects of turbulence (e.g. diffusive and non-diffusive 
effects) without the need to resolve the small scale turbulence. 
Non-diffusive effects include the a effect in mean-field dy- 
namo theory (Krause & Radler 1989) and the A effect in the 
theory of stellar differential rotation (Riidiger 1989). In order 
to calculate turbulent transport coefficients one usually makes 
the assumption that the equations for the fluctuating quanti- 
ties can be linearized (Moffatt 1978, Krause & Radler 1980). 
This implies that one is forced to restrict oneself to the case 
where the fluctuations are weak compared with the mean field. 
This is clearly not very useful for astrophysical applications. 
Nevertheless, mean-field theory is usually applied in the as- 
trophysically interesting parameter regime (e.g., Riidiger & 
Kitchatinov 1993, Kitchatinov et al. 1994), which should in 
principle be well beyond the regime of validity of such a lin- 
ear approximation. How is it then possible that the results are 
actually quite reasonable? 

The approximation under which higher order terms are ne- 
glected is known as the quasilinear, the first order smoothing, 
or the second order correlation approximation. With the excep- 
tion on an important additional term, the results derived under 
this approximation turn out to be similar to those derived un- 
der a more general approximation, which is sometimes referred 
to as the minimal tau approximation (MTA). This approxima- 
tion has its roots in early papers by Vainshtein & Kitchatinov 
(1983), Kleeorin et al. (1990, 1996), and its general usefulness 
has recently been stressed further with the papers by Blackman 
& Field (2002) and Radler et al. (2003). A recent review of non- 



linear dynamo theory and MTA has been given by Brandenburg 
& Subramanian (2005a). MTA has also been applied to the case 
of passive scalar diffusion (Blackman & Field 2003), where its 
validity has been tested numerically (Brandenburg et al. 2004). 

In the hydromagnetic case, the single most important ad- 
dition to the theory is the attenuation of the a effect by a 
term proportional to the current helicity of the small scale 
field. The beauty of MTA is that this term emerges in a com- 
pletely natural fashion. In the first order smoothing approach 
this term is absent, although even there it can be incorporated 
in a more phenomenological and hence less convincing fash- 
ion (Brandenburg & Subramanian 2005a). The role of the cur- 
rent helicity term was first identified by Pouquet et al. (1976) 
and later associated with catastrophic quenching by Gruzinov 
& Diamond (1994) and many others after them. The purpose 
of the present paper is to test MTA in the context of dynamo 
theory. 

2. The method 

In mean field dynamo theory one splits the magnetic field B 
into a mean magnetic field B and a small scale field b — B — 
B, and derives the mean-field dynamo equation 

— = Vx(UxB + £- v J). (1) 

Here £ = u x b is the turbulent electromotive force, J = 
V x B/hq the mean current density, /i the vacuum perme- 
ability (assumed unity throughout the rest of the paper), and 
the velocity U = U + u has also been split into mean U and 
small scale u = U — U velocities. The crucial difference be- 
tween the first order smoothing approximation and MTA is that, 
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instead of writing an expression for the electromotive force 
£ = u x 6, one considers the expression for its time deriva- 
tive, 



d£ -. =- 1 

— = uxb + uxb. 

at 



(2) 



The evolution equations for u and 6 contain terms that are lin- 
ear in the fluctuations. These terms lead to quadratic correla- 
tions of velocity (which are considered to be known from a 
turbulence model) and quadratic correlations of the magnetic 
field (which have to be obtained from a dynamical feedback 
model). In addition, there are triple moments resulting from 
the quadratic nonlinearities in the equations for u and 6. The 
sum of all triple moments is denoted by T and is assumed to 
be equal to the normalized quadratic correlations —£/t, i.e. 



— £ 

T = (closure hypothesis). 

r 



(3) 



Blackman & Field (2002) motivate this term on physical 
grounds: in the absence of any helical "driving", i.e. when 
the quadratic correlations vanish, £ should always decay. It is 
mainly the numerical simulations that can lend some support 
to this otherwise rather ad hoc and ill-justified assumption. [In 
principle one can also do this closure in Fourier space and adopt 
a k -dependent r(fc), as has been done by Radler et al. (2003) 
for example; this is discussed in more detail in Brandenburg & 
Subramanian (2005a).] 

As was first shown by Blackman & Field (2002), the use 
of MTA leads to an explicitly time-dependent equation for the 
electromotive force which, in the isotropic approximation, is 



d£ — - £ 

-r-=aB -r) t J 

at t 



(MTA). 



(4) 



This has to be contrasted with the usual result from the first 
order smoothing approximation (FOSA), 



£ = aB-r] t J (FOSA), 



(5) 



which agrees with Eq. (4) if the time derivative of £ is ne- 
glected. Here, a — to. is the usual a effect and i] t — Tij t is 
the turbulent magnetic diffusivity. One should point out, how- 
ever, that the meaning of t is different in the two cases, and that 
FOSA can only be justified if the r in this formalism is small. 

In the passive scalar case a relatively robust method for es- 
timating t from turbulence simulations is to impose a mean 
concentration gradient across the box, and to treat the devia- 
tions from this mean gradient as periodic in all three directions 
(Brandenburg et al. 2004). In the present case this corresponds 
to imposing a mean magnetic field, Bo = constant. The mag- 
netic field is then given by B = B + b, where b is the depar- 
ture from the mean field, b. Since dB /dt = 0, the evolution 
of b is given by 



6 = V X [U X (B + 6) - T]j + / mag ], 



(6) 



where a dot denotes a time derivative, j = V x b is the cur- 
rent density, and the current density is measured in units were 
/io = 1- We have allowed for the possibility of a magnetic forc- 
ing term, / mag , which is set to zero unless stated otherwise. 



We assume that the gas is weakly compressible and isothermal 
with constant sound speed c s . The evolution equation for the 
velocity is therefore 



u = -u ■ Vu - Vh + j x (J5 + b)/p + / kin + F v 



(7) 



where h is proportional to the enthalpy. Since the ratio of spe- 
cific heats is unity, we have h = h + c s In p, where ho is 
a constant whose value is unimportant for the dynamics. The 
evolution equation of In p is 



d In p j dt = —u ■ V In p — V • u. 



(8) 



We define a generic forcing function / (see Appendix A) 
and put /kin = / in the kinetically driven case (/ mag = 0), or 
/mag = / with /kin = in the magnetically driven case. The 
wavevector k(t) of the forcing function is delta-correlated in 
time, so at each timestep a new vector is chosen randomly with 
kf — | < \k\ < kf + i. We consider two cases: kf = 1.5 and 
kf = 5. The forcing function is chosen to be maximally helical 
with positive helicity. 

In the main part of this paper we adopt combined volume 
and time averages, so, B = Bq = const. This means that in 
Eq. (4) both J and d£ / dt vanish. Therefore we simply have 



tolBo- 



(9) 



Expressing £ in terms of a = (£ ■ Bo)t/B , where {...}t de- 
notes a combined volume and time average, we have simply 
a = ra. Here, a only depends on the quadratic correlations. 
Using the evolution equations (6) and (7), this leads to 



a 



I (-LJ ■ u + p Q 1 j ■ b) 



(10) 



This important relation was first obtained by Pouquet et al. 
(1976). In the following we will use the abbreviations 

a = a K + a M , &K = ~\u ■ u, a M = \Po 1 j ■ b. (11) 

The triple moment, T, has two separate contributions from the 
two evolution equations, T — Tm + Tk, where 



Tm = u x V x (u x 6), 



T K = {-u ■ Vu -Vh + jx b/p) x b, 



(12) 



(13) 



and 6 = 6 + B (1 — p/ Pa) includes the fluctuating contribution 
from the applied field due to density fluctuations. 

So far, and throughout most of this paper, we consider the 
case of an imposed constant mean field, B , with zero mean 
current, Jq = V x Bq = 0. While this is reasonable for mea- 
suring the value of a, it is of course unrealistic for dynamo 
calculations, where there must be a mean current proportional 
to the curl of the mean field. Furthermore, if a mean field is im- 
posed, this mean field can obviously not change by the dynamo 
action. Nevertheless, it provides a useful method for calculating 
a. This method gives similar results than the more complicated 
and less accurate methods that can be invoked in dynamo sim- 
ulations (Brandenburg 2001, see his Figs 14 and 15). However, 
in order to show that Eq. (3) has any justification, we insert a 
discussion (Sect. 3) where we consider the local correlation in 
a simulation without an imposed field. 
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Fig. 1. Comparison of the spatial dependence of two components of 
the mean magnetic field and the triple correlation in Run 3 of BO 1. The 
magnetic field is normalized by the equipartition field strength, _B eq , 
and the triple correlation is normalized by fciB^q, but scaled by a fac- 
tor of 2.5 to make it have a similar amplitude as the mean field. Note 
that B x (solid line) correlates with T x (filled dots) and B z (dashed 
line) correlates with T z (open dots). 

3. Spatial variation of T 

Before we present in detail results for the case of an imposed 
field, we consider first the case of a dynamo-generated mag- 
netic field. In that case the magnetic field shows a marked si- 
nusoidal variation in one of three possible directions. Which of 
the three directions is chosen depends on chance. In the fol- 
lowing we compare with a simulation of homogeneous heli- 
cally forced turbulence of Brandenburg (2001); in his Run 3 
the magnetic field shows a slow variation in the y direction, so 
we define mean fields by averaging over the x and z directions. 

Figure 1 shows that the x and z components of T have 
a positive correlation with the mean magnetic field. Note that 
£ itself has a negative correlation with B, for a negative a. 
And since this simulation has positive kinetic helicity, the a 
obtained is indeed negative. Therefore the positive correlation 
of T with B is consistent with Eq. (3) with a positive r. (The 
term t] t J does not change the sign of aB, but only reduces its 
magnitude. This follows from the fact that in a dynamo sim- 
ulation £ has to overcome the magnetic diffusion term.) We 
have used B as a proxy for £ in order to study the correlation 
of T with £, because £ itself is even more noisy than T(y). 
This is because aB almost cancels (77 + %) J and the fact that 
J consists of derivatives of B which contributes further to the 
noise. 

In conclusion, our first assessment of T is very encouraging 
in that its phase relation with the magnetic field is such that it 
would indeed contribute as a damping term, i.e. that it has the 
expected sign. 

4. Reynolds number dependence of a and a 

Having shown that T varies locally in the expected sense, we 
now concentrate on the case of a uniform imposed field where 



it is possible to use averages over the full box. (In the dynamo 
case considered in the previous section such averages would 
have given zero.) We solve Eqs (6)-(8) numerically. The forc- 
ing amplitude is such that we are in the low Mach number 
regime, i.e. |«| <C c s , so for all practical purposes the flows can 
be considered nearly incompressible. We consider runs with 
different hydrodynamic and hydromagnetic Reynolds numbers, 

Re = u ims /{vkf), R ia = u ims /(rik{), (14) 

but keep the magnetic Prandtl number unity in all cases, i.e. 
R m = Re. Here, kf is the wavenumber of the energy-carrying 
scale which is assumed to be the forcing wavenumber of ei- 
ther 1.5 or 5. The simulations have been carried out using the 
Pencil Code 1 which is a high-order finite-difference code 
(sixth order in space and third order in time) for solving the 
compressible hydromagnetic equations. In the following we 
consider four different values of R m and increase the resolution 
correspondingly by factors of 2 between 64 3 and 512 3 mesh- 
points, and in some cases between 32 3 and 256 3 meshpoints. 

We measure velocity in units of the sound speed c s , density 
in units of the average density p , magnetic field in units of 
yfjlopoCs, and length in units of fc-f , where k\ is the small- 
est possible wavenumber in the box. Whenever possible we 
present the result in normalized (explicitly nondimensional) 
form. 

We begin by considering the dependence of a on R m ; see 
Fig. 2. We focus on the results for two different field strengths, 
Bq = 0.03, corresponding to Bo/B cq ss 0.45, and Bq = 0.1, 
corresponding to B n /B cq w 1. Throughout this paper, error 
bars are obtained by taking the extrema from 5 separate aver- 
ages, each over 1/5 of the full time series. 

In both kinetically and magnetically driven cases a is neg- 
ative when the helicity of the forcing is positive. It turns out 
that, when the turbulence is kinetically forced, —a decreases 
with increasing R m like where n is between 1/2 and 1. 
For large field strengths and large values of R m we expect an 
asymptotic dependence with n = 1 (Cattaneo & Hughes 1996), 
but when these values are not yet large enough the data are best 
fitted with values of n less than 1 . 

In the magnetically driven case a is generally, at compa- 
rable values of w rms , larger than in the kinetically driven case. 
Furthermore, in the magnetically driven case a does not show 
the systematic decrease of R m seen in the kinetically driven 
case. This can be understood as a result of a modified Keinigs 
(1983) relation; see Appendix B. 

Next, we consider the values of &k and ccm in both ki- 
netically and magnetically driven cases; see Figs 3 and 4. The 
results depend on the values of Bo and fcf. In the kinetically 
driven case there is a clear tendency for ctM to increase with 
i? m until it approaches and slightly exceeds the value of — oik- 
For larger values of R m , —5k and «m are constant and close 
to each other. 

In the magnetically driven case, oik and au depend only 
weakly on i? m . Furthermore, qm can exceed — <Sk by even an 
order of magnitude; see the first panel of Fig. 4. 

1 http : / /www . nordita . dk/ software /pencil- code 
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cases where the turbulence is driven kinetically with a body force /kin 
(lower graph) or magnetically with a forcing term / m a g in the induc- 
tion equation (upper graph). The vertical bars on the data points give 
estimates of the error (see text). 

It is not easy to present the data of all runs in a meaningful 
way in a single graph. This is why we have summarized the 
results of all runs in tabular form; see Table 1, allowing alter- 
native representations and interpretations to be made. 

5. The relaxation time 

After these preparations we can now return to the main question 
addressed in this paper: what is the relevant relaxation time in 
MTA and how does it depend on R m . Assuming a steady state, 
we can use Eq. (4) and calculate 



r = a/(&K + «m)- 



(15) 



The value of r obtained in this way is shown in Fig. 5, where 
we show the results separately for the kinetically and magnet- 
ically driven cases. The general trend is quite similar to the 
dependence of a on R m . 

As we have discussed above, at large R m , &m can even 
slightly exceeds the value of — &k- Since a itself does not 



change sign, this would imply that r would become nega- 
tive for large R m . An obvious solution to this problem is that 
(5k and &m should be preceded by some additional quenching 
functions, gx and gu, respectively. Both functions are expected 
to be of the order of unity, but they may not be exactly equal; 
see Kleeorin et al. (2002). In the following we combine the in- 
formation from the kinetically and magnetically driven runs to 
compute separately t<7k and t<?m, and assume qk = rg^a^ 
and aM — rgM&M- We thus have an equation with two un- 
knowns, 



a = Tg K a K + rg M aM, 



(16) 



where the two unknowns are rg K and rg M . 

In order to obtain a second equation, we can consider the 
simulations for kinetically and magnetically driven turbulence 
as two independent measurement. This leads to two separate 
measurements of a, &k, and q.m, distinguished by superscripts 
(k) and (m), respectively. In this way, Eq. (16) becomes a ma- 
trix equation 



*( k ) 
„("») 



~(k) 



~(k) 



K 



( TgK 
\Tgm 



(17) 
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Table 1. Summary of the measured rms velocity and normalized as for kinetically and magnetically forced runs, together with the resulting 
Strouhal numbers and their corresponding errors obtained by taking the extrema from 5 separate averages over 1/5 of the full time series. 
Values in parenthesis show departures from the trend and should be regarded as uncertain. For kf = 1.5 the resolution varies between 64 3 and 
512 3 meshpoints for r\ = 2 x 1CT 3 and 2 x 10~ 4 , respectively, while for kf = 5 the resolution varies between 32 3 and 256 3 meshpoints for 
rj = 5 x 1CT 3 and 5 x 1CT 4 , respectively. 
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When the field strength is increased, both tcxk and t<5m 
are quenched uniformly, i.e. approximately independently of 
R m . This becomes clear when comparing the first two panels 
of Fig. 6, or by inspecting the data of Table 1 . For a fixed value 
of R m (around 50) the quenching with increasing field strength 
is shown in Fig. 7. This quenching, which sets in once B /B cq 
is of the order of unity, is rather strong (proportional to Bq 3 ). 
Similarly strong quenching has previously been found for a if 
the contribution from j ■ b is ignored (Moffatt 1972, Riidiger 
1974, Riidiger & Kitchatinov 1993). 

6. Energy and helicity spectra 

In Figs 8 and 9 we present, both for the kinetically and mag- 
netically driven cases, shell-integrated spectra of magnetic and 
kinetic energies, M(k) and E(k), respectively, as well as cur- 
rent and kinetic helicities, C(k) and F(k), respectively. They 
are normalized such that / M(k)dk = \{B 2 ), j E{k)dk = 
\{u 2 ), J C(k)dk = (J ■ B), and / F(k)dk = (u ■ u). 
The two helicity spectra are subject to a realizability condi- 
tion, \C(k)\ < 2kM(k) and \F(k)\ < 2kE(k) (Moffatt 1978). 
If the bound is nearly saturated, we say that the field is fully 
helical. It turns out that for k > fc f we have \C(k) \ < 2kM(k) 
and <C 2kE(k), so magnetic and velocity fields are in 

fact not fully helical at small scales. 

Since both C{k) and F(k) seem to scale like fc~ 5 / 3 in 
the inertial range, both helicities are governed by contributions 
from large scales. This implies in particular that the small scale 



which can be solved for t^k and t<?m- However, since the 
root mean square velocities, n£ms and u£mL are different 
in the kinetically and magnetically driven cases, we non- 
dimensionalize each measurement independently, so we define 

(k,m) = a( k, m)/u (k, r ) 5 (lg) 



U K,M — a K,M / ^f" 



(k,m) 
rms 



(19) 



St.g K ,M = MrmsfcfTgK,M, (20) 

where u rms = [w^msUms] 1 / 2 is the geometrical mean of the 
rms velocities for kinetically and magnetically driven runs. 
Thus, the matrix equation that we actually solve is 



« (m V V4 m) i m) / v St 5M 



(21) 



The result is shown in Fig. 6. Note that, unlike the previous 
case, t^k and rgM are now always positive. In some cases the 
results are quite similar to the previously determined values of 
r, but in the case with _B = 0.03 and k f = 1.5, where r 
became negative (not seen in the logarithmic representation in 
Fig. 5), both t<7k and rgM are now positive, and approximately 
constant for R m < 100. However, for larger values of R m , both 
values may decrease, although it should be noted that the error 
bars are also larger and the accuracy of the error itself may not 
be reliable either. Thus, it is not yet clear that the decline seen 
for the largest value of R m is indeed real. 
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Fig. 4. Dependence of a^ 1 ' and a^' on _R m in the magnetically 
forced case. Vertical bars give error estimates. 

current helicity, (j ■ b), and hence also 5m, should not depend 
on the length of the inertial range and hence should be indepen- 
dent of R m . If this were not the case, i.e. if the magnetic field 
were fully helical with a fc~ 2 / 3 current helicity spectrum, (Sk 
would be asymptotically proportional to rU 4 (Blackman & 
Brandenburg 2002). Our results suggest that this is not the case 
here. This seems to be in conflict with earlier results for helical 
dynamos (Brandenburg 2001), but those results were at 4 times 
smaller resolution (120 3 instead of 512 3 meshpoints) and the 
forcing was at 3 times smaller scale, so the effective length of 
the inertial range was 12 times smaller and hence completely 
absent. What remains true, however, is that the field is nearly 
fully helical at the energy carrying scale, i.e. at wavenumber 
k { . 

7. Kinetically vs magnetically forced turbulence 

In order to gain more insight into the nature of the turbulence 
in kinetically and magnetically forced cases we show visualiza- 
tions of the field-aligned components of velocity and magnetic 
field, u z and B z , respectively. In the kinetically driven case we 
see a large scale pattern in velocity, consistent with power at 



Fig. 5. Strouhal numbers as a function of R m for different combina- 
tions of Bo and fcf. Here, gn = gm is assumed. The horizontal lines 
are drawn to indicate the range over which the Strouhal numbers are 
approximately constant. In the second panel therefore, the slope of the 
lower graph is —0.68. 



the expected energy-carrying scale corresponding to the forc- 
ing wavenumber fcf = 1.5; see Fig. 10. The magnetic field, 
on the other hand, is dominated by smaller scale structures. 
The strength of the imposed field is still less than the equiparti- 
tion value (Bo/B cq w 0.4) and there is no obvious anisotropy. 
Although the velocity field is helical, there is no obvious large- 
scale pattern in the magnetic field. There are two possible rea- 
sons for this: lack of scale separation (see Sect. V of Haugen 
et al. 2004) and suppression of large scale dynamo action by 
the large scale field (Montgomery et al. 2002, Brandenburg & 
Matthaeus 2004). 

The magnetically driven case is quite different. There 
seems to be a large scale pattern in the magnetic field (Fig. 1 1). 
Although there is no scale separation, it is possible to drive 
a larger scale field in this case. The main difference is that 
in the kinetically driven case, dynamo-generated large scale 
and small scale fields have opposite sign of current helicity 
(Brandenburg 2001). In the absence of sufficient scale sepa- 
ration this leads to cancelation. In the magnetically driven case 
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e -2 / 3 , where e is the total energy dissipation rate. The two helicity 
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all scales, they would be close to the corresponding energy spectra, 
which is not the case. Instead, both energy and helicity spectra show a 
short subrange where they are best fitted with a k~ 5 /' A power law. The 
helicity spectra, compensated by k 5 ^ 3 , are shown in the lower panel. 
The turbulence is kinetically forced, with imposed field Bo = 0.1, 
using 512 3 meshpoints. 
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the current helicity has the same sign at all scales, so there is 
no cancellation and therefore the build-up of a large scale field 
is possible. However, when the field is too strong (Bo = 0.1, 
corresponding to Bo/B cq ps 1) no large scale pattern in the 
magnetic field develops; see Fig. 12. 

8. Conclusions 

Although a closure hypothesis of the form T = —S/t 
has a plausible physical interpretation as a relaxation term, 
(Blackman & Field 2002), it is nevertheless quite crude and 
lacks a rigorous physical basis. Nevertheless, on empirical 
grounds this closure assumption is quite appealing. First of all, 
unlike the first order smoothing approximation, MTA provides 
a natural and convincing explanation for the j ■ b correction 
term to the a effect in Eq. (10). This is indeed a key ingredi- 
ent to the dynamical quenching model (Kleeorin & Ruzmaikin 
1982), where the evolution of j ■ b is obtained by solving the 
magnetic helicity equation for the small scale field. The dynam- 
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Fig. 9. Same as Fig. 8, but for the magnetically forced case with Bo = 
0.1, 512 3 meshpoints. 



ical evolution of the j ■ b term provides a good model for the re- 
sistively limited saturation (Field & Blackman 2002, Blackman 
& Brandenburg 2002, Subramanian 2002), which is seen in 
helically forced simulations in periodic boxes (Brandenburg 
2001, Mininni et al. 2005). Another appealing property of MTA 
is that the additional time derivative, d£ / dt in Eq. (4), restores 
causality and prevents infinitely fast signal propagation. In the 
case of turbulent passive scalar diffusion, one can easily see 
that this term turns the heat equation into a damped wave equa- 
tion, where the maximum signal speed is the rms velocity of 
the turbulence (Brandenburg et al. 2004). 

One of the striking features of the dynamical quenching 
model is that the resistively limited saturation and the resis- 
tive quenching of the effective a can already be explained 
without invoking any explicit R m dependence of r. Instead, 
the i? m dependence of the steady state values of a emerges 
solely as a result of magnetic helicity conservation and might 
hence be alleviated if there is a helicity flux that offsets mag- 
netic helicity conservation (Blackman & Field 2000, Kleeorin 
et al. 2000, 2002, Vishniac & Cho 2001, Subramanian & 
Brandenburg 2004, Brandenburg & Sandin 2004, Brandenburg 
2005, Brandenburg & Subramanian 2005b). However, there is 
no compulsory reason that there might not be an explicit i? m 
dependence of r after all. Our present simulations begin to shed 
some light on this question. However, the results are still not as 
clear cut as one would like them to be. Looking at Fig. 6, it is 
evident that for R m < 100, St is independent of R m and de- 
pends only on Bq\ cf. the uniform decrease of St between the 



first and second panels of Fig. 6. However, for R m > 100, St 
may actually decrease with R m , although the error bars also in- 
crease. If this decrease is indeed real, it might indicate that the 
closure hypothesis adopted here is too simplistic. 

We should emphasize that, although we have in this work 
focused on the a effect, this is not the only effect that produces 
large scale dynamo action. An example is the shear-current ef- 
fect of Rogachevskii & Kleeorin (2003, 2004). Here the elec- 
tromotive force has a component in the direction of — W x J, 
where W = V x U is the vorticity of the mean flow. (In the 
present work this effect is of course absent because there is no 
mean flow and our mean fields are defined as volume averages, 
so J = 0.) Our work is relevant even in this nonhelical case be- 
cause, like the a effect, the W x J effect produces mean fields 
that are helical. Depending on the amount of magnetic helicity 
losses from open boundaries, there will be a tendency to limit 
the production of net magnetic helicity within the domain, so 
small scale current helicity of opposite sign must be produced 
at the same time as large scale field is generated. The presence 
of small scale current helicity implies that there must be a fi- 
nite ayi- This is an example where there can be an cum term 
even though «k = 0. Related examples where qk = 0, but 
om 7^ 0, are the reversed field pinch (e.g. Ho et al. 1989), and 
the so-called selective decay (e.g. Montgomery et al. 1978), as 
described in Blackman & Brandenburg (2002) and Yousef et 
al. (2003). Another example where — 0, but «m ^ 0, is the 
Vishniac & Cho (2001) mechanism whereby dynamo action is 
accomplished through the current helicity flux that drives chm\ 
see Brandenburg & Subramanian (2005b) for model calcula- 
tions showing that this effect produces saturation field strengths 
independent of the magnetic Reynolds number when the flux 
exceeds a certain threshold. 

In the present paper we have only discussed the case where 
the helicity of the turbulence is driven explicitly. Although this 
is a common assumption that allows progress to be made in 
that the flow can be taken to be statistically isotropic, it should 
be noted that in astrophysical applications helicity is usually 
introduced by the interaction between stratification and rotation 
(Ki-ause & Radler 1980). To calculate this effect, MTA has to 
be used once more (e.g. Radler et al. 2003), which then renders 
the a effect proportional to r 2 . Nevertheless, as is emphasized 
in Brandenburg & Subramanian (2005a), also in this case the 
a effect is attenuated by the current helicity contribution from 
the small scale field. 

Earlier simulations of helical dynamos (no imposed field) 
suggested that the field is bihelical, i.e. with one sign of helic- 
ity at large scales, corresponding to wavenumber k\, and the 
opposite sign at small scales, corresponding to wavenumber kf 
(Brandenburg 2001). If the field were to continue to be nearly 
fully helical even at smaller scales, this would have implied 
that q;m would be attenuated by a RU 4 factor (Blackman & 
Brandenburg 2002). Our present results suggest that this is not 
the case. In agreement with both closure calculations (Andre 
& Lesieur 1977) and direct numerical simulations (Borue & 
Orszag 1997) we find approximate /c~ 5 / 3 scaling for the ki- 
netic helicity (see also Ditlevsen & Giuliani 2001). Again, 
this implies that the velocity field is not fully helical beyond 
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Fig. 10. Visualizations of u z (left) and b z (right) in the kinetically forced case. Bo = 0.03, v = n = 2 X 10~ 4 , kf = 1.5, 512 3 meshpoints. 




Fig. 11. Visualizations of u z (left) and b z (right) in the magnetically forced case. Bo = 0.03, v = T] = 2 x 10~ 4 , k t = 1.5, 512 3 meshpoints. 




Fig. 12. Same as Fig. 11, but for Bo = 0.1. There is no clear evidence for a large scale field. 
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wavenumber k t ■, so a K should also be independent of i? m . This 
is indeed the case, except for very small values of R m . 

There are various ways of extending the present studies of 
turbulent transport to other fields. An obvious possibility is tur- 
bulent viscosity and the A effect for modeling the Reynolds 
stress tensor (Kapyla et al. 2005). For example, one would 
expect the relaxation time to depend on the angular velocity 
(Kitchatinov & Riidiger 1993). This dependence could be mea- 
sured in a similar way as the dependence on the magnetic field 
studied in the present paper. Another related application is the 
study of passive scalar transport in the present of rotation. Both 
aspects are quite essential for many astrophysical applications. 
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Appendix A: The forcing function 

For completeness we specify here the forcing function used in the 
present paper 2 . It is defined as 

f{x, t) = Re{Nf k{t) exp[i*(t) ■ x + i#t)]}, (Al) 

where x is the position vector. The wavevector k(t) and the random 
phase — 7r < <f>(t) < ty change at every time step, so f(x,t) is 6- 
correlated in time. For the time-integrated forcing function to be inde- 
pendent of the length of the time step St, the normalization factor N 
has to be proportional to St^ 1 ^ 2 . On dimensional grounds it is chosen 
to be N = f c s (\k\c s /St) 1/2 , where f is a nondimensional forc- 
ing amplitude. The value of the coefficient fa is chosen such that the 
maximum Mach number stays below about 0.5; in practice this means 
fo = 0.01 . . . 0.05, depending on the average forcing wavenumber. At 
each timestep we select randomly one of many possible wavevectors 
in a certain range around a given forcing wavenumber. The average 
wavenumber is referred to as k{. Two different wavenumber intervals 
are considered: 1...2 for kf = 1.5 and 4. 5. ..5. 5 for fcf = 5. We force 
the system with transverse helical waves, 

/fe = R-4 n ° hel) with R ij = 5 Ji^0*L, (A.2) 

where a = 1 for positive helicity of the forcing function, 

/<r hol) = (fe x e) /Vfc 2 -(fe-e) 2 , (A.3) 

is a non-helical forcing function, and e is an arbitrary unit vector not 
aligned with k; note that |/j^| 2 = 1. 

2 This forcing function was also used by Brandenburg (2001), but in 
his Eq. (5) the factor 2 in the denominator should have been replaced 
by y/2 for a proper normalization. 



Appendix B: Modified Keinigs relation 

The R m dependence of a can generally be understood as a direct 
consequence of magnetic helicity conservation (Gruzinov & Diamond 
1994). Dotting Eq. (6) with the vector potential a, where b = V x a, 
and adding the corresponding evolution of a ■ b, averaging over a pe- 
riodic volume, and using (u x Bo) x b = — (u x b) ■ Bo we obtain 

i^(a ■ b) = -{u x b) ■ B - n(j -b) + {b- / mag ). (B.l) 

Using (u x b) = aBo, we find in the steady state 

a = - v ^ + ^pi. (B.2) 

The first term, which was already obtained in an early paper by 
Keinigs (1983), is manifestly ^-dependent [see also Brandenburg & 
Matthaeus (2004) for a derivation of this term], but the second term is 
not. This explains the behavior seen in Fig. 2. 
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